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SUMMARY 

Equations are derived for computing the form of an oblique shock 
■wave as it passes through supersonic regions in which static pressure, 
stagnation pressure, stagnation temperature, or combinations of these 
are continuously variable. Rigorous portions of the analysis are limited 
to shock strengths for which the flow downstream of the shock remains 
supersonic. , When no downstream waves other than those generated by 
the interaction process are present, the rate of change of shock angle ’ 
with upstream Mach number is found to he a function only of the local 
shock angle and upstream Mach number ^ hence, the prapagatlon through 
a nonuniform region depends only on the initial shock strength and Mach 
number. A procedure is described for computing the supersonic portion 
of the flow field downstream of the shock wave. 

For the special cases of supersonic shear flow and Prandtl-Meyer 
flow, charts of the shock eingle as a function of upstream Jfech number 
are presented so that the passage of a shock wave throu^ these types 
of nonuniform regions can he easily traced. For a prescidhed initial 
shock strength and, initial Mach number, a minimum upstream Mach number 
is foiond below which no physically realistic solution can be obtained 
with the equations for simple propagation. This result serves as a 
sufficient condition for the avoidance of separated flow, reversed 
flow, or other upstream effects. An example is computed of the propaga- 
tion of a shock wave through a wake-type supersonic shear profile and 
the flow field downstream of the shock is constructed. 


INTRODUCTION 

The most frequently discussed example of the propagation of a 
shock wave into a nonuniform flow region appears, at present, to be 
the shock-wave boundary- layer interaction problem. Although much 
theoretical and experimental research has been devoted to this problem, 
it is still unsolved in the sense that no method exists whereby the 
effect of a shock of prescribed strength on a boundary layer with 
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prescribed profile can be predicted. The difficulty in formulating a 
theoretical approach can be attributed chiefly to the fact that the 
flow upstream of the shock is not independent of the presence of the 
shock. If the subsonic portion of the boundary layer -irere removed 
(by continToouB suction, for example), it might be expected that this 
difficulty would be eliminated. Circumstances frequently arise, how- 
ever, for which even this eaqjectation is not fulfilled. It was shown 
qualitatively in reference 1, for example, that the presence of a 
supersonic trake upstream of a blunt body requires modification of the 
upstream flow in order that physically realistic pressure gradients 
can be obtained downstream of the detached shock. The fact that the 
flow upstream of the shock is completely supersonic does not, therefore, 
guarantee that the upstream flow will be uninfluenced by the presence 
of the shock. 

In order to formulate criteria for deteamiinlng whether or not the 
flow upstream of a shock is independent of the presence of the shock, 
it is first necessary to analyze the propagation of shock waves throu^ 
regions in which upstream flow conditions are variable over a \ride 
range. Such an analysis was carried out for weak waves and for the 
particular case of a supersonic shear layer in reference 2. 

In the present report, the propagation of oblique shock waves of 
arbitrary strength through supersonic regions in which static pressure, 
stagnation pressure, stagnation temperature, or combinations of these 
are variable over wide ranges is considered. In addition to providing 
criteria for the avoidance of upstream effects of shock ^mves, the 
a na l y sis provides a method for tracing shock -vraves throixgh jets or 'vrakes 
of .know profile or throtigh other types of flow with large nonunlforml- 
ties. An analysis of the propagation of weak waves in general nonuniform 
regions is also given, and a procedure is discussed for constructing 
the supersonic portion of the flow field downstream of shock waves or 
weak waves. This investigation was conducted at the HACA Lewis 
laboratory. 


ANALYSIS 

General. Equations 

The interaction of a shock wave with a nonuniform supersonic stream 
is first discussed for the most geneial conditions feasible. The result- 
ing equations are then simplified for several special types of flow. 
Symbols used are defined in the appendix. 

For the most general case considered, the following assumptions 
are made ; 

(l) In the vicinity of the interaction region, effects of viscosity 
turbulence are negligible. This assimiption permits use of nonviscous 
shock and flow equations. 
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(2) The flow is everywhere supersonic and, consequently, changes 
in pressure and flow direction taJse place only across shock ira-ves or 
characteristic lines (idiach waves) . 

(3) The flow field can he divided into stream tubes of small thick- 
ness. Stagnation temperature T and ratio of specific heats y are 
constant in each stream tube but may change across the streamlines that 
boimd each stream tube, (if the static t^perature, and consequently 

y, changes appreciably across the shock wave, an average value of y 
can be used, for each stream tiibe with slight additional complication. ) 

(4) Stagnation press\ire remains constant in each stream tube except 
when that stream tube passes through the shock wave. 

(5) Fluid properties upstream of the shock -vrave contaln.no large 
discontinuities, (if this condition is not ‘imposed, the numerical work 
soon becomes prohibitively lengthy. Evaluation of the exact solution 
at the point of impingement of a shock wave on a contact surface that 
separates fluids whose properties differ by large amounts is discussed 
in references 3 and 4.) 

With the preceding assumptions, the t/ave and streamline pattern 
near a point on the shock wave in the nonuniform flow field can be 
represented as in sketch 1. Lines Oa and Oc are Incident waves or 
characteristics upstream of the shock; line Og is an incident wave or 
characteristic in the downstream field; and Of is a reflected wave or 
characteristic. Static pressure emd flow direction may change across 

any characteristic or across 
the shock wave; stagnation 
pressure can change across 
any streamline or across 
the shock wave; and stagna- 
tion temperature and speci- 
fic heat can change across 
any streamline. Assump- 
tions (3} and (5) imply that 
all waves except the shock 
wave are sufficiently weak 
so that the relation bebvreen 
flow deflection and pressure 
change across the wave satis- 
fies the weak- wave equation: 

where primes represent down- 
stream values and unprimed 
symbols denote upstream 
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values. The plus and minus signs axe used, respectively, for 'VTaves 
whose slope is positive or negative in sketch 1. All waves except the 
shock iTave are inclined at the local Mach angle relative to the local 
flow direction. 

The following relations "between values upstream of a shock (unprimed) 
and values downstream of the shock (primed) are required for the analysis 
(see reference 5) : 

— = fi(M,r:.<P) = — ^ M2sin2<J) - lUL-i ‘ (2) 

p i^'"'’/y- + i' r + 1 


X = f2(M,T,<p) tan 


-1 


O 

M sin q> cos q> - cot (p 
1 + _ sin^q)^ 


f3(M,T;.<P) = 


_ yM 


,2 


M’^ - 1 


where 


(t + l)%^sin^ - 4(l^sin^q> - l) (YM^sin^<P + l) 
M = 

[gyM^sin^ (p - (y - i) J [[(y - l)M^sin^T + 2^ 


C3) 

(4) 

(5) 


The equations relating the quantities in 
regions of sketch 1 are derived ■^d.th reference 



each of the several 
to the numtering shorn 
in sketch 2, which is 
an expanded version of 
the vicinity of point 0 
in sketch 1. The waves 
and the streamline 
separate the flow near 
point 0 into eight 
regions, in each of 
which all fluid proper- 
ties are constant. Con- 
ditions in the upstream 
regions 1, 2, 3, and 8 
are kno-vmj and the inci- 
dent shock angle (pg is 

known. From the shock 
equations, conditions in 
region 7 are determined 
■by (pQ, Mq, and yg. 

The unknOTOS are cp^ 


Sketch 2 
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and conditions in regions 4, 5, and 6. Since pressure and flow angle 
are constant across streamlines, the following relations exist among 
the pressures and the flow angles: 


^4 “ ^5 ( 6 ) 

( 

•P3 Pg P7 Pg P3 ^ 

Substitution from equations (l), ( 2 ), (3), ( 4 }, and (s) into equation (?) 
yields 




1 - 






C^i - \)\ 




But from equation (3) 


^4 “ ^3 ^ 2 ^%^'>^ 3 ’>‘’^ 3 ) 


(9) 


X 7 a Xq + 

Hence, equation (8) becomes 


(lOj 




raMs^ 


L 


(Xj - Xg + f2(%'T5»<P5) 


|i + ^3(“a'Tra'‘%5[^6 ■ ^a ■ ^zC“a'Te''*’a]J^i(^''>'a'’’a) ^ 


All quantities in equation (u) are known except ip^. Mg, and Xg. if 

I'raves of type Og in sketch 2 can be feliminated from the problem, the 
solution is' greatly simplified, since then X^ a X^ and Mg ^ M^j and 

the only remaining unknown is the transmitted- shock angle <p 3 . Condi- 
tions under which waves of type Og can be neglected can be ascertained 
with the aid of the following sketches; ' 
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In sketch 3(a) waves of type Og sire generated "by the object that 
produces the shock wave and reach the shock before it passes through 
the nonuniform, flow region. For this case^ therefore, this type of 
incident wave is obviously not negligible. In sketch 3(b), on the other 
hand, waves of type Og that originate on, or are reflected from, the 
object reach the shock only after it has passed through the nonuniform 
region. There remains the possibility, however, (as illustrated in 
sketch 3(b)) that waves reflected within the nonuniform field itself 
■^d.11 approach and appreciably influence the shock form before the shock 
leaves the nonuniform flow region. Conditions for which these waves 
are negligible can be determined by an analysis of the strength of the 
^mYe of type Og that separates regions 2 and 3 in sketch 3(b) relative 
to the strength of the wave of type Of that separates regions 1 and 2. 
Thus, in sketch 3(b), since p^ = Pj, p^ » Pj^, X 3 => and =: Xg, 

use of equation (l) yields f with — J^L__ _ f (m) ) 

\ Vm^ - 1 J 

CX 3 - Xg) [f(M 2 ) + f(l%) - f(Mg) f(Mi)(X 2 - Xi)]=x (Xg - X^) [f(Mi) - f(Mg)] 
from which is obtained 

V 

^ -^2 _ f(Mp) - f(%) ^ ^ 

f(M 2 ) + !(%) - f(Mj^) f(M 2 )(X 3 - X^^) “ (2 - f^X)f 
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Consequently, ireives of type Og, if they arise from internal reflection 
in the region of nonuni form flow, can he neglected without serious error 
if fAX « 2, where AX' can he Interpreted as the strength of the strong- 
est I'rave of type Of. Since the strength of the Of waves is determined 
hy the magnitude of the discontinuities upstream of the shock wave, this 
criterion, in effect, specifies the magnitude of the upstream discontinu- 
ities that can he tolerated if equation (ll) is to he simplified. It 
I'/ill he assumed throughout the remainder of the report that the upstream 
flow satisfies this requirement, and that consequently M0 and V pnu 

he set equal to My and Xy, respectively, if no waves other than those 

resulting from shock pi-opagation are present. The modifications reqxiired 
when craves of type Og other than those resulting from internal reflec- 
tion are present are considered in the section DISCISSION AND AETLICATION. 


If Xg =i X^ and Mg « My, equations (S) or (ll) become 










[H 









?3 


( 12 ) 


Since changes in flow variahles from one stream tube to the next 
have been assumed small, functions of Mg, CP3, and pg differ only 

slightly from the same functions of Mg, <Pq^ and pg. Consequently, let 


^n “ 1^2,3) 


and let 


Pg =■ Pg + dp 

^3 “ ^ 


Equation (l2) then becomes 


fl + df 1 


P8 

Pq + <ip 


[1 - fgCdX + dfg)] 


(13) 


,(14) 


“ 1 - f ^ 2 'Q Cl5) 

where the differentials correspond to changes in quantities upstream 
of the shock \7ave from one stream tube to the next in the direction of 
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propagation of tlie shock wave. If only first-order differential terms 
are retained, equation (l5) becomes 


df-i 


= - f3(dX+ dfg) - ^ 


(16) 


But 


dfn =* 

Sfl 

d(p + i 



dy 

X 

S<P 


br 

dM 

df 2 = 

afg 

d<P + j 

r^f2 

Sfp 

+ 

dr 

S<p 

^Sm 

hr 

dM 


dX 

dX = ^ dM 


dp 

dp = dM 
dM 


) 




dM 


dM 




(17) 




Hence, equation (l6) becomes 

3f 


d<p 

dM 


1 

f 


L. ^ 7 

1 ^ dM/ ^ \ 


2m 




dM 


t)-- 

dM/ p 


1 ^^1 ^^2 
fjL ^<p 


B tp 


dM 


(18) 


All quantities in the right member of equation (l8) depend only 
on conditions upstream of the shock ira.ve; hence, if the variations of 
y, X, and p -;ri.th upstream Mach number are specified (as they v/ill be 
if upstream conditions are knovjn), then the rate of change of shock 
angle viith upstream Mach number is a function only of (p and M. For 
any initial shock angle, and for arbitrary (but continuous) upstream 
Mach number distribution, eqixatlon (l8) can therefore be integrated to 
obtain the form of the shock ira-ve in regions of nonuniform supersonic 
flow. For discontinuous, but small, changes in upstream conditions, 
equation (l8) must be evaluated in a step-by-step manner from one 
discontinuity to the next. The validity of equation (l8) is limited 
to shock waves for which the flow downstream is supersonic and to 
segments of the shock vrave that are free of interaction vrLth downstream 
incident waves other than those resulting from the propagation process. 

The partial derivatives required in equation (18) as evalioated from 
equations (2) and (3) are: 
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AyiM sln^(p 


2YM^sin^(p - (y - l) 


1 ^^1 4yM^sin «p cos <p 


2rM^sin2q> - (y - l) 




2(M^sin^q> - l) 


Sm 


^^2 
S (p 


3f2 

hr 


^ (y + l) [2YM^sin^(p - (y - l)] 

(y + 1 )m cot <p 


|l + ^ ^ - sin2<pj + cot^<p (M^sin^cp - l)^ 


- ^ -- | ^ - 2 sin2<p^ ^ - r sln^cp^ M^sin^(p 


sin' 


■<p'^ 1 + ^ ^ - sin^cp^J + cot^<p (M^sin^(p - l)^ 


- ^ cot (p (M^sin^<p - l) 


"* T 3 

1 + ^X_i_i _ sin^tp^ + cot^9 (M^sin^tp - l)^ 


(19) 


( 20 ) 


( 21 ) 


( 22 ) 


(23) 


(24) 


Special Cases 

Supersonic shear flov . - If flo-w angle^ static pressure, and ratio 
of specific heats are constant immediately upstream of the shock, equa- 
tion (is) reduces to ^ ^ 

1 ^fl hf^ 

dq, f^ InT ^ 

-^=f4(M,9) (25)- 


(3M 


. _= 
* ^5 5 ,,, 


A contour plot of dcp/dM as function of cp and M is sho^jn in 
figure 1. Integral ciirves of figure 1 are presented in figure 2. The 
construction of these integral curves hy an isocline method is aided 
hy the fact that the ciirve of Mach angle against M (M sin (p = i) is 
itself an integral curve of equation (25) . This fact can he deduced 
from equation (25), which reduces for this case to 
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'dqA _ tan <P 1 

sin<P =1 “ ■ M “ " 


(26) 


wMcli is exactly the slope of the curve M sin <P = 1. Consequently, 
the slope of all integral curves at their intersection with an 
f^ = constant contour is the same as the slope of the M sin <p •= 1 
curve at its intersection with that contour. 


The integral cu 2 rves of figure 2 yield the variation of shock angle 
■vjith upstream Mach numher when the reflected wave from each differential 
interface is taken into account. For comparison, figure 2 also contains 
curves obtained i-Tith the assumption that these reflected waves are 
negligible and, hence, that the pressure ratio across the shock (or 
M sin <p (see equation (2)) remains constant as the shock passes through 
the shear region. Since the slope of any M sin <p a constant curve is 
d<p/dM =» - (tan <p) /M, the curve f 4 (M, <p) = - (tan <P)/m (fig. 2) is the 

locus of points for which the slope of the integral curves is equal to 
the slope of the constant-pressure-ratio curves. Comparison of the 
constant-press-ure-ratio cirrves with the integral ciurves in figure 2 
shows that the pressure ratio across the shock decreases as upstream 
Mach number decreases in the region to the left of the curve 
f 4 (M,T) ts - (tan <P)/m. To the ri^t of this curve, pressure ratio 
increases as upstream Mach nxnnber decreases. The intersection of this 
curve •t'llth the Mach angle curve at M = 1.41 agrees with a res\ilt 
obtained in reference 2 for the interaction of weak waves \)lth parallel 
supersonic shear layers. For a weak incident wave entering the region 
of decreasing Mach number, it ras foimd that for M » -y/2 no reflected 
rave occured. The curve f 4 (M,‘P) = - (tan <p)/m in figure 2 thuB rep- 
resents a generalization of this result to the case of raves of arbitrary 
strength. 


For fairly weah shocks, the assumption that the pressure ratio 
across the shock remains constant through the shear region appears to 
be an adequate approximation. For stronger shocks, however, the pres- 
sure gradient downstream of the shock can become quite large, and the 
actual shock angles may differ considerably from those obtained rath 
the constant-pressure-ratio assumption. 

The integral curves of figure 2 cannot be constructed beyond the' 

M' » 1 curve, since the reflected-rave concept used to derive these 
curves is meaningless if the dorastream flow is subsonic. For this 
region (cp^ < cp<90°), however, it should be approximately valid in 
certain cases to assume that the pressure Immediately downstream of 
the shock is constant along the shock, particularly if the flow down- 
stream of the shock is expected to be almost parallel. A large pressure 
gradient normal to the subsonic streamlines would be inconsistent with 
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tlie exj)ectecL nearly paoi’eillel flow. Conaequently, in order to compute 
the form of the shock for cp > (p^, the integral curves of figure 2 can 
probably be extended to cp ts 90° along curves of constant-pressure 
ratio with no significant error. 


■When such extensions are made, it is seen ftom figure 2 that for 
any initial Mach number and shock angle there exists a minlmnm upstream 
Mach nimber below whicih no shock solution is obtained. If the Madi 
number profile of the shear region is such that the ttiI nlnmni Mach number 
is less than that for which a shock of ' specified initial strength can 
be traced through the region, some readjus-tment of the upstream flow 
appears to be required to ob-tain a solution that is physically realistic. 
For example, the situation sho'wn in sketch 4 may be considered, wherein 

‘ a shock of initial angle 

cpQ 28.6° and initial 
Mach nuniier of 3.0 enters 
a shear region in which 
the miniTiTum Mach number 
is 1.2. From figure 2, the 
incident shock is found to 
become normal at M = 1.36. 
If the shock remained nor- 
mal from a to b, as 
indicated by the dotted 
extension of the shock, 
then the pressure ratio at 
b would be 1.52 as com- 
pared -with the pressure 
ratio at a of 2.0. Since 
the normal shock implies 
initially parallel subsonic 
do-wnstream flow, it seems 
unreasonable to expect that 
a pressure gradient normal 
■to the streamlines such as 
required in sketch 4 will 
be established. 



A situation similar to that sho'wn in sketch 4 -was discussed in 
reference 1, where the modification required in the shock form when a 
shear layer occurs ahead of a blunt body -was considered. It -was argued 
that the pressure near the base or center line of the shear layer could 
be increased to a magnitude compatible trith that outside the shear 
layer only if the stagnation point moved upstream to the vicinity of 
the shock -wave. Hence, a separated-flow region -was required ahead Of 
the body to pro'vlde a mechanism for establishing pressure equidJLbrium 
in subsonic portions of the do'wnstream flow. Although the nature of 
the modification required in the present case is less evident, the 
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formation of a separated- flow region iri-tli a cx>nseq,uent larntda shock 
outside of the shear layer appears to be the most likely possibility. 

In Einy event, the curves of figure 2 provide a criterion for determining 
whether simple propagation of a shock wave throu^ the shear region 
can take 'place. If the shock can be traced throi:i^ the shear region 
by means of figure 2, then the physical shock should correspond to that 
calculated. If, on the other hand, the shock cannot be traced through 
the shear layer, then a more complicated flow pattern will be obtained 
which may or may not involve reversed flow and flow separation.®’ This 
criterion yields no Information, of course, on the magni’fcude of the 
separation that may result if simple propagation is impossible. 

Honunlform Isentropic flow . - If the stagnation pressure upstream 
of the shock is constant, then from the isentropic flow relation 

p = P ^ “T"" 

there ia obtained 


i ^ _ yM (1 + 

p dM 


m2)- 


(27) 


If X is also cons-fcant, equation (l8) becomes 


dtp 

dM 


A 

f 


+ ^3 ^ (l + m2) 


1 ^ ^■'-2 

St * 


5!f 

St 


(28) 


t- 

N 

lO 

N 


Since the ■variation of X with M depends, in general, on the structure 
of the flo-tf upstream of the shock, the profile of M and X must, in 
general, be specified -to obtain the shock form from equation (28) . For 
some types of flow, however, such as Prandtl-Meyer flow, the quantity 
dl/dM is a prescribed function of M so that equation (28) can be 
represented on a T - M plot as for shear flow. Thus, for Prandtl- 
Meyer flow, since 


®’Since spearated-flow regions unsupported by solid bomdaries are 
not commonly observed experlmen-fcally, it may be helpful to point out 
that such regions can be maintained in equilibrium relative to the main 
stream by the high dotmstream pressure that is established in the sepa- 
rated region. The branched- shock configuration that results from flow 
separation can produce the required equilibrium conditions. 
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* Vr?i 


(M^ 1) + sin-1 I + c 


there is obta±ned 


<3M 


± -^/l 


- 1 


M ^1 + 


(29) 


and eq.uation (28) "becomes 

Sf. 


dp 

dM 


1 ^ . "''2 
f 1 ^3 'Sir 


1 

Stp ^ 


3 3(p 



= fsCM.^P) 
(30) 


where the plus sign is used when dX/dM is positive and the minus 
sign when dX/dM is negative. 

Positive values of d^dM correspond to interactions such as 
those shown in sketches 5(a) and 5(h); negative values occur for inter- 
actions such as those shown in sketches 5(c) and 5(d) and arise less 
frequently in practice than positive values. 

The compressive, or reversed, Prandtl-Meyer flow shown in 
sketches 5(h) and 5(d) is, of course, described hy equation (29) only 
in regions in which the flow can he represented hy a sequence of weak 
waves, t/hen the weak waves form an envelope shock of sufficient 
strength to produce notable entropy increases, the flow angle and Mach- 
number are no longer related hy the Prandtl-Meyer equation. 

The first term of the right member of equation (30) is the value 
of d(p/dM for constant p, X, and y (equation (25)); the second 
term is the effect of' static pressure and flow an^e changes in Prandtl- 
Meyer flow. 


A contour plot of dcp/dM = f5(M,<p) for Prandtl-Meyer floi-r 

is shown in figure 3 for the case in which the shock is entering a 
region in which both X and M either increase or decrease 
(dX/dM >0). The curve sin<p = 1 /m is again an integral curve of 
equation (30), since the last term In this equation is zero when <p is 
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(c) Decreasing increasing M (d) Increasing \ decreasing M 



Sketch 5. - Types of interaction of shock "waves "with Prandtl-Meyer 

flow regions. 

the tiach angle. Comparison of figures 1 and 3 shows that for most 
■values of M and <p the shock angle increases less rapidly with 
decreasing Mach number for Prandtl-Meyer flow than for shear flow. 

Integral curves of figure 3 are plotted in figure 4, "together 
"vfith several constant-pressure-ratio curves. It is seen that the 
assumption "that pressure ratio across the shock remains cons"tant in 
passing thro"Ugh the nonuniform region is even less "valid for this case 
than for shear flow. For fairly weak shocks (X* _ X < 10°), however, 
it is found that the integral curves follow quite closely ciorves of 
constant- flow deflection (cons"bant fg)* The integral curves for 

shear flow, on the other hand, followed curves of constant-pressure 
ratio more closely than ciirves of constant-flow deflection. For very 
weak; "waves, "the "two "types of curves hecome indistinguishable. 
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The curve f 5 (M,< 5 >) *= - (tan <P)/m separates the region in which 

pressure ratib decreases iTith decreasing Mach number from the region 
in which pressure ratio increases as Mach number is reduced.. This 
curve again intersects the Mach angle curve at M = -y/2. 


Strength Distribution for Weak Waves 


As mentioned previously, an analysis of the interaction of weak ■ 
waves with a supersonic shear flow was developed in reference 2. In 
this analysis, the assumptions were made (as in the preceding analysis 
of shock interaction) that the flow upstream of the wave was uninflu- 
enced hy the presence of the wave and that this upstream flow could 
be divided into elementary stream tubes. An additiohal assimiption 
was made that the Mach number in each stream tube was unchanged in 
passing through the wave. In the present section the interaction of 
a weak wave 'td.th more general nonuniform supersonic regions is consid- 
ered, and the assumption that the Mach number remains constant in each 
stream tube is omitted. 


A weak wave is herein defined 
and flow deflection are related by 



as one across which pressure ratio 
equation (l) < 1 -Jhen the upstream 
floif variables are specified, the 
local downstream flow variables 
are determined by means of a five- 
sector configuration (sketch 6 ) 
in which, since P 3 = p^, the 

pressures are related by 

^ ( 31 ) 

P2 Pi P5 Pi 

Use of equation (l) for the pres- 
sure ratios of equation (3l) 
results in 


1 + 


r ±^2 
■7^2^ “ 


“^ 2 ) 


!-• 


ri^i 


L - 1 




P5 

^1 


1 -- 


- 


(^4-^5) 


(32) 


As for the shock wave, an average value of y instead of the upstream 
value can be used for each wave with slight eidditional complication. 

In this mariner the variation of y with temperature can be taken into 
account as the computation proceeds. 
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Since a.i T clianges across the -waves and across the streamlines are 
small, let 

P5 = Pj_ + 4Pj_ H p + dp 
f (n^Mg) = + Sf(ri.Mi) = f + &f 

f(r5^J%) = + if a f + ±f 


■where 


Equation (32) can then he -wri-tten (since ^4 = 

|1 + (f + 8 fK>i . + Xg - Xi + Xi - X2)] [1 - f(X2 - \)] - (1 + §) [1 - (f + af)(x^ - Xg)] (33) 

Let dX» X5 - and let 9 he the change in flow angle across the 
weak -wave; then 

Q ts Xg - Xj^ 

d0 = (X^ - X^) - (Xg - Xj^) 


f(r,M) = 


tm^ 


Vm^ - 


0 + d.0 — 


Equation (33) nov ‘becomes 

jl + (f + Sf)(e + d0 + dX - 0)j jl - foj = - (f + df)(0 + d0^ 


which, when rearranged, hecomes 


f aBja + ^^i + y^-fe^i + y) + Y + ^]--eeLf + ^(i-ef-ecLc)- ax(i - ef)(f + sf) 
Neglect of an terms of higher order than the first yields 
f d0(2 - 0f) - 0 if + - f dX^l - 0f) 


( 34 ) 
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Calculations show that f is less than 10.1 for 1.01 < M < 5.0. 
Consequently, in this range of Mach numbers, if 6 is of the order 
of 1°, then 6f can he considered negligible compared with unity and 
eqixation (34) becomes 


e 2f 20f Vp j 


For the case of shear flow (dp =0, dy = 0, 
reduces to 


e 


- 2 dM^ 

4(M^ - 1). 


(35) 


d\ = O) , equation (35) 


(36) 


which is the result derived by a different process ih reference 2 with 
the assumption that Mach nimber is constant in each stream tube. As 
pointed out in reference 2, it is readily shown from equation (36) that 
a weak incident wave entering a region of decreasing Mach number pro- 
duces a reflected wave of the same sense (compression or expansion) or 
oprosite sense depending on whether M is less than or greater than 
->/2. Integration of equation (36) yields, as in reference 2, 


A = ^ f - '1 

^ \Mq 2 - ij 


(37) 


where subscript zero refers to any initial or specified values. The 

fxinction (M^ - 1)^^^/M is plotted in figure 5. For shear regions 
with large changes in Mach number, it is evident that the cumulative 
effect of the reflected waves, which accounts for the change in 0 in 
the interaction process, is negligible only if the initial weak wave 
is itself negligible. 

For problems in which dp/p or dX are of the same order as 
0 df , equation (35) miist be used instead of equation (37) . If the flow 
upstream of the weak wave is of the ,Ik-andtl-Meyer type, eqixation (35) 
becomes, by use of eqmtions (27) and (29), 

M^ - 2 tM (l ± l) 

e 2 m(m 2 - 1) ■ + l_|J=M2y • 20f 

where the plus sign corresponds to dX/dM > 0 (see sketches 5(a) 
and 5(b)). Since an incident weak wave is everywhere inclined at the 
local Ifech angle, interactions of the type shown in sketches 5(c) 
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and 5(d) do not occur in tMs case (d6 ^ 0 'because dM *= O) . For 
the physically significant case (dX/dM > O) equation (38) can he 
written 


1 ^ 
9 dM 


- 2 

2M(M^ - l) 


+ 


0M 


- 1 

(l + 


(38a) 


This equation is easily integrated numerically for any prescribed 0 q 
and 1^. 


DISCUSSION AND AmJCAPION 


Construction of downstream flow by stepwise process . - Men the 
form of a shock waYe In a nonuniform region has been determined by 
eqmtion (l8) or its corollaries, and when the resulting downstream 
flow is completely supersonic, the method of characteristics for anisen- 
tropic flow can be used to complete the -computation of the downstream 
flow field. This method, however, is q.ulte tedious numerically; and it 
may be desirable, in cases where great accuracy is not required, to 
establish the essential features of the downstream flow by a simple 
method. The weak-wave procedure to be described constitutes such a 
simplification. 

Before discussing the computation procedure, it is necessary to 
establish and define the order of magnitude of the various waves 
encountered. For this purpose, the interaction process shown in sketch 7 
is considered, wherein a shock 'irave is shoim entering a nonuniform 
region from the upper left. The flow upstream of the shock may be 
continuous or it may contain discontinuities sufficiently small so that 
the criterion for neglect of all except the primary reflected waves for 
computing the shock form is satisfied. This criterion was discussed in 
the section General Equations. After the shock form is. determined, the 
upstream flow field can be divided into a sequence of streamtubes in each 
of which all fluid quantities can be assumed to be constant. This sub- 
division is required for stepwise computation of the downstream flow 
field even if the upstream flow is continviDUB. The subdivision produces 
a ifave pattern of the type shown in sketch 7. 
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The shock wave can he considered to he a zero-order wave which, in 
passing throng the nonuniform region, produces a sequence of reflected 
first-order waves. The integrated effect of the first-order waves on 
the strength of the shock wave can he of the same order of magnitude as 
the initial strength of the shock. A zero-order solution to the inter- 
action piHDhlem must therefore retain the integrated, effect of first- 
order waves. Similarly, each of the reflected first-order waves, whose 
initial strength is known from conditions just downstream of the shock 
■vrave, can he considered to he a primary weak wave passing through a 
region of nonuniform flow. Each first-order wave produces a sequence 
of second-order reflected waves, whose integrated effect on the strength 
of the first-order waves can he of the first order. ’ To compute accu- 
rately the strength of the first-order waves as they emerge from the 
noniiniform flow field, it is therefore necessary to use equation (35) 
for each first-order wave (since p, X, and M are all variable 
upstream Of these waves) . The use of equation (35) to compute the 
strength of each first-order wave yields an approximate solution for 
conditions upstream of the next first-order wave. With these conditions, 
the strength distrlhutlon of each succeeding first-order wave can he 
computed. 
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Sketcli 8 stiDVB tlie precise manner in "wtiich tlie second-order waves 

are taken into account "by this pro- 
cedure. Conditions upstream of first- 
or^r reflected waves ah and fd 
are known from the shock form. The 
strength of the secondary reflected 
wave he is obtained from equa- 
tion (35), which gives the change 
in flow deflection of irave ahh ' 
at h. The strength of wave he, 
in turn, determ in es conditions 

upstream of segment de, from which 

the strength of segments de and dg 
can he calculated hy equation (35) . 

By continuation of this procedure, 
the entire supersonic portion of the 
downstream flow field can he con- 
structed. The only approximation 
made is that waves such as he do 
not significantly affect the average 
flow variahles upstream of primary imve segments such as fd; in other 

ac 



words, that the ratio 


is close to unity. Since this ratio is 


ad 

very close to unity when sill -vraves are inclined at the local Mach angle, 
the maximum error that can result from violations of this condition is 
equivalent to the error that could result from neglect of a single 
sequence of second-order waves (those nearest the shock ivave), which, 
in turn, is approximately equivalent to the effect of a single first- 
order -t-Tave. In general, therefore, the use of equation (I8) and its 
co3Xjllaries for shock waves and the use of equation (35) and its corol- 
laries for first-order waves will yield an accurate first-order solution 
to the interaction problem. 


In many problems, the foregoing procedure may yield higher accuracy 
than is reqiilred; and it may he desirable, for the sake of simplicity, 
to Ignore second-order waves entirely. In such a case, the pressure 
ratio or flow deflection across each first-order wave is considered 
constant until it intersects another first-order wave or a shock wave. 
The neglect of second-order '^•raves is permissible, if only a fenr such 
waves are generated. Th'us, in sketch 7, the first three or four first- 
order waves generate only a few second-order waves before they emerge 
from the nonuniform region. In general, if the upstream Mach number 
veiries over a small range, so that only a few streamsheets are required 
to break the flow into stiff! ciently small differences, then the neglect 
of second-order waves will result in negligible error. 


Effect of incident first-order traves . - If first-order waves other 
then those arising from the passage of a shock t/ave through the nonunifomn 
region are present, then the computation of the downstream flow field by 
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the wave method Is feasible only if second-order waves are entirely 
neglected. Even with this appi-oximation^ however, the precise construc- 
tion of the flow is rather laborious because the intersections of waves 
with waves, and waves with streamlines occur at different positions. 

The situation locally is as shown in sketch 9, where a, b, and c may 
be first-order waves such as those shown in sketch 3(a}j d, e, and f 
may be first-order reflected waves due to the passage of a shock wave; 

and g and h are stream- 
lines. The simplest computa- 
tion procedure consists in 
assianing that stagnation 
pressure and stagnation 
temperatures are consta nt 
in grids such as j kmn 
(rather than between stream- 
lines) . If the width of 
each grid is of the same 
order as the initially 
assumed streamtube width, 
then no appreciable error is 
introduced by this procedure. 
The streamlines can be 
ignored except insofar as 
they indicate the approximate 
locad stagnation pressure 
and temperature, and each 
computation of properties in 
successive grid spaces is 
merely an application of the weaE-wave relation (equation (l)) to the 
problem of interaction of two weak waves. Thus, if conditions in grid 
spaces 1, 2, and 4 are known from previous computations, conditions in 
region 3 are obtained from the relation 



Sketch 9 


P 5 P 4 ^ P 5 Pg 
. P 4 ^’l ^2 '^1 

or, with the use of equation (l) for these pressure ratios. 


(39) 


[1 - f(M 4 )(X 3 - \J]^ + f(%)(X 4 - Xi)] - [1 + f{M 2 )(X 3 - Xgj] [l - f(lli)(X 2 - X^]] 
which, when solved for kg, yields 

~ - X^)] + f({'^^)X^[l + f(Mj^)(X^ - X]_)] + f(M3^) [Xg - Xj^ + - Xj 

f(l%] [1 + - XJ] + fCMg) [1 - f(Uj^)(X 2 - x^)] 
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Since all -vraves are consl^red to "be weak, let 

fCMg) =» f(Mi) + Bf ‘ 
f (M^} =« f (M^^) + df 

If the steps are sufficiently small so that f(M)AX is small compared 
with unity, then neglect of all except the lowest order terms in equa- 
tion ( 40 ) yields 


— kg ^ ^4 "* ^1 
or 

kg - k^ a kg - k^ 

which shows that the assumption commonly made for isentropic flow, 
namely, that the flow deflection remains constant across each first- 
order wave in interactions with other first-order waves, is also appli- 
cable when continuoiis or stepwise small variations in entropy and energy 
occur. 

With P 3 and 13 taken equal to P]_ and x±> respectively, 

equations (4l) permit calculation of pressure and Mach number in 
region 3. The accuracy of this procedure depends on the negligibility 
of second-order reflections at the differential entropy and energy 
interfaces -which, in t\im, depends on the number of stream tubes required 
to divide the upstream flow into sufficiently small differences. If the 
upstream Mach number varies over a very -wide range, uae of the method 
of characteristics for anisentropic flow may be required in place of 
the first-order wave procedure if a high degree of accuracy is required 
downstream of the shock wave. 

Interaction of incident first-order waves with shock wave . - When 
first-order waves other than reflected waves due to passage of the shock 
are present, it is likely that some of these waves will intersect the 
■ shock wave in the region of nonuniform flow. This possibility was 
pointed out in the first section of the AHALYSES, hnd discuaslon was 
postponed because the analysis neglecting the waves iTas valid if such 
intersections were few and for segments of the shock between such inter- 
sections. The local modification of shock strength required at these 
intersection points is derived with the aid of sketch 10 , in which ein 
incident wave of type Og is assumed to intersect the shock wave 
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■between two successive refiected 
waves due to upstream flow vari- 
ations. At the intersection- 
pointy an additional reflected 
wave Of can arise which is 
due only to the interaction of 
Og with the shoci; wave and not 
to variation in the upstream 
flow. The strength of Og at 
the , Intersection point is con- 
sidered known from computations 
of its interactions with all 
first-order waves in its path. 
Although a solution for the 
strength of the wave Of result- 
ing from the interaction with 
the shock; wave can he obtained 
from equation (ll)^ the procedure 
is rather labortous, and it seems 
adequate for most purposes to 
Sketch 10 neglect the reflected wave Of 

entirely. The assunq)tion that 
no reflected wave of type Of 

occurs at the intersection is the same assumption that is made with good 
accuracy in- the shock- eacpansion theory for two-dimensional airfoils. 

The conditions for which the Of wave is negligible compared with the 
incident wave are. fully discussed in reference 6 . Specifically, the 
assumption is that in sketch 10 Xg = and pg « P 3 , where pj and 
Xj are known from the preceding step in which the formula for the 
intersection of two weak waves was used. The change in pressure ratio 
across the shock at point 0 transfers the shock to another integral 
curve of equation (iS) (or fig. 2 or 4), which it follows Tintll another 
Intersection of this ,ty^ occurs. The entire shock form is therefore 
determined to first order by its own propagation process and by first- 
order waves from external sources that Intersect it. 

The accuracy with which the flow field can be calculated by neg- 
lecting all except first-order wves and the initial shock wave depends 
on the strength of the second-order waves that are neglected. An 
estimate of the error can be made by determining the strength of the 
second-order waves that are produced by the strongest first-order wave 
in the flow field. The errors in the computed pressures and flow angles 
should be comparable with the changes in pressure and flow angle across 
these neglected waves. 
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EXAMPLE 

As an example of shock propa^tion •through supersonic shear regions^ 
the form, of a shock -wave in a constant-pressure ■wake-'type profile 'tra.s 
computed. The profile assumed is sho-wn in figure 6 , -together -with the 
shock form "traced -withi the aid of figure 2. The initial shock Eingle 
and Machi number -were assumed -to he 34.0° and 3.0^ respectively. The 
streamlines and first-order reflected -waves -were computed under the 
assumption that second-order -waves were negligible and that consequently 
the flow deflection and pressure ratio across each first-order reflected 
wave remain cons"tant. 

The TniT Tl muiii Mach number in the -wake -was -taken to be 1.75, which 
yields a region of subsonic flow downstream of the shock, but does not 
exceed the criterion for simple propagation. The region of subsonic 
flow was desired "to illus-trate an approximate procedure for handling 
such regions when they are imbedded in a supersonic stream. The pro- 
cedure used was as follows: In the region of subsonic flow, the wave- 

type calculation^ of course, is not "valldj but it -was assumed that the 
pressirre difference and flow deflection across each first-order wave 
was impressed on the imbedded subsonic region and that these pressure 
differences propagate through the subsonic region normal -to the flow 
direction and emerge at the opposite sonic line in the foim of waves 
with the same strength as the incident waves. Obis process appears 
reasonable if the extension of the shock through the region q> > (p^ by 
means of the constant-press-ure criterion is approximately valid. In 
that case, the shock form determines the first-order -wave pattern which, 
in "turn, determines the streamwise pressure gradient in the supersonic 
streamlines adjacent "to the subsonic streamlines. The assumption that 
this pressure gradient is transmitted normally thro-ugh the subsonic 
region is consistent -with the expec-fcation that curva-ture of the subsonic 
stream-lines will be small if the curvature of the adjacent supersonic 
streamlines is small. 

With the assumption that no separation takes place downstream of 
the shock and that pressure is ti’ansmitted normally through the subsonic 
portion of the dD-wnstream shear layer, the interaction process is easily 
constructed in terms of the local Mach angles and flow angles. The 
resulting pattern of streamlines and first-order reflected waves is 
shown in figure 6 for the example chosen. The flow angle just downstream 
of the shock -raried from 16.7° at the point of incidence to 19.2° near 
the center line. The pressure ratio across the shock varied from 3.12 
at the point of incidence -to 2.72 neeir the center line. Since the 
static pressure and flow angle are constant between any -tiro reflected 
first-order waves (due to neglect of second-order -waves), these figures 
represent the range of -variation for the entire interaction region. 
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Although these -rarlatlons in flow angle and pressiire ratio arfe 
appreciable relative to the initial flow deflection and pressiare ratio 
(about 15 and 13 percent, respectively), they ane nevertheless fairly 
small and may be negligible for some purposes. The resulting variations 
in Mach number are indicated on figure 6. In the stream outside the 
shear layer the variation is only about 4 percent of the initial down- 
stream Ifech number. These ranges of variation can, of course, be pre- 
dicted as soon as the form of the shock is known. The location of the 
variations, however, must be determined by constructing the first-order 
TOve pattern. 

The self-consistency of the method used to trace changes in pres- 
sure and flow angle through the subsonic streamtubes can be checked to 
some extent by comparing the change in streamtube area required by one- 
dimensional equations with the change actually obtained. In the present 
case, the change in Mach number from 0.90 to 0.80 along the central 
streamtubes requires an area change of only about 3 percent, which is 
self-consistent with the area changes actually obtained within the error 
in tracing the streamlines. 


COnCLUDIUG REMABES 

Analysis of the propagation of shock waves and first-order waves 
through nonuniform regions of supersonic flow showed that the strength 
distribution of an incident wave in such regions is determined by the 
upstream distribution of Mach number, flow angle, pressure, and specific- 
heats ratio provided that no traves other than the incident wave and 
reflected waves generated by its propagation are present. If the 
upstream flow can be divided into a sequence of small discontinuities, 
each incident wave can be considered to produce a sequence of reflected 
waves whose strength is an order of magnitude less than that of the 
primary wave. The integrated effect of these reflected waves, however, 
may be of the same order of magnitude as the strength of the primary 
wave. 


Men the primary wave is a shock wave, a minimum upstream Mach 
number is found below which no physically realistic solution can be 
obtained with the equations for simple propagation. A sufficient con- 
dition is thereby obtained for the avoidance of flow separation, 
reversed flow, or other nonsimple phenomena. 

Men first-order waves other than those reflected within the non- 
uniform region due to passage of a shock wave occur downstream of the 
shock, a first-order solution to the interaction problem is obtained 
by neglecting all second-order reflected waves . The accuracy of this 
procedure depends on the magnitude of the variations in fluid properties 
upstream of the shock and becomes less accurate when these variations 
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are large. If a high, degree of accuracy is required, tlie flow down- 
stream of tlie shock can he constructed hy the method of characteristics 
for anisentropic supersonic flow. 

Althou^ the wave analysis used becomes invalid when a portion of 
the flow downstream of the shock wave becomes subsonic, approximate 
procedures are given for computing the downstream flow field, as well 
as the shock form, when the curvature of the streanJjLnes in the down- 
stream flow field is small. 


Lewis Flight Propulsion Laboratory 

ITational Advisory Committee for Aeronautics 
Cleveland, Ohio, March 13, 1952 
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APPENDIX - SIMBOLS 


The foUowing symbols are used in this report; 


slW - 1 

fpCM^r^q)) pressure ratio across shock, p'/p (equation (2)) 
f‘ 2 (M^r^‘p) flow deflection across shock, - X) (equation (3)) 


fsCM^r.tp) 


YM 


,2 


■yjw 


(equation (4) ) 


- 1 


d<p/dM for shear flow (equation (25)) 

<3fl»/dM for Prandtl-Meyer flow (equation (30 ) ) 
Mach number 
stagnation pressure 
static pressure 
ratio of specific heats 


M 

P 

P 

r 

e 


<p 

<Pa 


flow deflection aqross shock or weak wave (X'' - X) (posi- 
tive value for compression, negative value, for expansion) 

flow angle relative to reference direction (positive sense 
indicated where needed) ' 

shock angle relative to local upstream flow direction 

value of (p for M' = 1.0 


For X, p, and M unprimed values represent conditions upstream 
of shock or weak wave and primed values refer to conditions downstream 
of shock or weak wave. 

Number subscripts for X, M refer to values in cor- 

responding regions in accompanying sketches. 
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Shook angle, t, deg 

Figure 2, - Variation of shook angle with upatroaii Maoh numher. Conatant upstream ntatlo preaaure and flow 

angle. Ratio of speolfio heats, 1.40., 
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Figure 3. - Rato oT ohango of nhook angle Kith upetnaan Mach nxupber In Prandtl-Hevor flow redon. Ratio 

of specific h«at8»*1.40. 
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